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DIFFUSION HITTING TIMES AND THE BELL-SHAPE 
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Abstract. Consider a generalized diffusion on R with speed measure m, in the natural scale. It 
is known that the conditional hitting times have a unimodal density function. We show that these 
hitting densities are bell-shaped if and only if m has infinitely many points of increase between the 
starting point and the hit point. This result can be viewed as a visual corollary to Yamazato’s 
general factorization for diffusion hitting times. 


1. Introduction and statement of the result 

This paper deals with a certain distributional property of hitting times for generalized diffusions. 
We use the standard notation for the latter processes, as described e.g. in Chapter V of [ 8 ]. Let m 
be a string, that is a right continuous non-decreasing function from [—00, +00] to [—00, +00], with 
m(—00) = —oo,m(+oo) = +00 and m( 0 —) = 0 . Set 

n = sup{x < 0 , m{x) = —00}, r2 = inf{x > 0 , m{x) = +00}, 
and define the positive measure m{dx) on [—00,+00] by 

m{dx) = dm{x) on (ri,r2), m([ri,r2]'^) = 0 and m({ri}) = m({r2}) = +00. 

Let {i?t, t > 0 } be a linear Brownian motion and {Lf, t > 0 , x € M} be its local time. Introducing 
the additive functional 

At = [ Lt m{dx) 

Jk 

and its right-continuous inverse Tt = inf{tt > 0 , Au > t}, we define the process 

Xt = Br„ 0 <t< C, 

with lifetime C = inf{t > 0 ,Xt = ri or r2}. Here and throughout, the notation inf 0 = -boo is 
implicitly assumed. The process X is strongly Markovian with state space Em = (xi, r2) fl Supp m, 
and we set for its law starting from x € Em- The measure m is called the speed measure of X. 
As is well-known, the above time-changed Brownian motion allows to construct all linear diffusions 
on an interval up to some monotonous, scale transformation - see again Chapter V. 7 in [ 8 ]. One 
can also choose a speed measure with discrete support, which leads to a class of continuous time 
Markov chains called gap diffusions in the literature - see [ 6 ] for details and examples. The above 
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generalized diffusion is chosen on the natural scale, which does not cause any loss of generality in 
our problem. 

The hitting time of y G Em is defined by 

Ty = inf{t > 0] Xt = y}. 

This dehnition extends to y = r* for i = 1,2 if we suppose \ri\ < oo and G Em, setting 
Tr^ = Ty In this situation we adjoin rj to Em and denote by Em the extended state space. 

Let X G Em and y 7 ^ x G Em be such that Px[tj/ < 00 ] > 0. Without loss of generality we will 
suppose y > X. This paper deals with the conditional hitting time distribution 

[Ty G dt] 


Txy[dt^ — 


[Ty < OOj 


It is well-known that the probability measure iTxy is absolutely continuous, and we denote by fx,y 
its density function on (0,-t-oo). It was shown in Theorem 1.2 of [S] that this density is always 
unimodal, in other words that it has a unique local maximum. 

In this paper, we are interested in the following refinement of unimodality for fx^y A smooth 
density function defined on a real interval is said to be bell-shaped if all its derivatives vanish at 
both ends of the interval and if its n—th derivative vanishes exactly n times, for all n > 1. Setting 
n = 1 shows that a bell-shaped density is strictly unimodal. For n = 2 the bell-shape property 
entails, as for the familiar bell curve, that there is one inflection point on each side of the mode 
and that the second derivative is successively positive, negative, and positive. The visual meaning 
of the bell-shape for n = 3 or 4 is less immediate and we refer to the introduction of m for details 
and references. Some standard density functions like the Cauchy, the Gaussian, the Gumbel or the 
Student are bell-shaped, as can be seen by a direct analysis. Showing this property for less explicit 
densities is however a more demanding task. 

Gonsider now the generalized inverse Gaussian distribution, with density 


g 3 ;A-lg-(xa; '-+'ipx) 


over ( 0 , -|-oo), where c is the normalization constant and the variation domain of (A, x, V’) is described 
in [1] p. 49. If X > 0, the above density is bell-shaped, as can again be seen by a direct analysis, 
and it is also a hitting density for some SDE driven by Brownian motion - see Theorem 2.1. in [T]. 
If X = 0, the above density is not bell-shaped, and it is the hitting density of a generalized diffusion 
only if A < I - see Gorollary 2 in |I2] . In this case it is also completely monotone, hence the hrst 
passage density to the nearest neighbour of some continuous time Markov chain with discrete state 
space - see [2] p. 147 and the references therein. If x > 0 the speed measure of the underlying 
diffusion has an everywhere positive density, whereas if x = 0, A < I, the involved speed measure 
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is atomic. This example suggests a general relationship between the bell-shape property for the 
hitting times of a generalized diffusion and the support of its speed measure. This connection is 
illustrated by the following characterization. 

Theorem. The density fx^y is bell-shaped iff m has an infinite number of increase points on {x,y). 

This simple criterion applies to all SDE’s driven by Brownian motion, and also to more singular 
processes like those of Examples 2.12 and 2.13 in [3]. There exists an immense literature on hitting 
times for real diffusions, however it seems that the above interesting distributional property has 
not been investigated as yet. Our method to prove the theorem relies on a general factorization 
by Yamazato [12], and a total positivity argument which the first author had used in a previous 
paper m on stable densities. In this respect, it is worth recalling that diffusion hitting times are 
always infinitely divisible. The bell-shape property for all positive self-decomposable densities with 
infinite spectral function at zero is stated as an open problem in im - see Conjecture 1 therein. 

2. Proof of the Theorem 

2.1. Proof of the if part. Translating the measure m if necessary, we may and will suppose x = 0 
and y > 0. We appeal to the factorization obtained in Theorem 1 of which reads 

T^Oy = hi * h2, 

where hi has Laplace transform 



for an at most countable family of positive parameters {oj} which is either empty or increasing and 
such that h 2 is an absolutely continuous probability with completely monotone 

density function. Observe that the probability distribution /ii belongs to the so-called Bondesson 
class - see Chapter 9 in m- The Levy measure of the infinitely divisible distribution p 2 also fulfils 
a certain condition which is described in the Theorem of m, but this will not be used in the sequel. 
Let us finally refer to Section 15.2 in uni for a thorough presentation of this classical result. 

From the proof of Theorem 1 in [T2| p. 155, we also know that 

r= Th) 

where f>{z^ A) is the unique continuous solution to 

(1) = z X {z — t)fi{t, \)dm{t), z > 0. 

A[0,2) 
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In particular, the family {—Oj} is the set of zeroes of the function A i-a A). We first show the 
intuitively obvious fact that under the assumption that m has infinitely many points of increase on 
(0,y), the family {oj} is inhnite. This simple remark is more or less already included in |10] p.204, 
but we write it down for completeness. 

Fix n > 1 and let 0 < ?/i < ... < Vn+i = y be such that each interval (y^, yi+i) contains a point 
of increase of m. For each i G [1, n] and A > 0 we have 

rvi+i 

i’ivi+iA) = i’iviA) + / 

> (yi+i - y*)-^ (y*, A) 

> A(yi+i-yi)i/)(yi,A)m(yi,yi+i) 

where the inequalities follow from ([T]) and the fact that z i-A ip{z,X) is positive, increasing and 
convex on (0, y). An induction argument shows that there exists a positive constant K such that 
V’(y) A) > KX^ for all A > 0. This being true for every n > 1, we see that A i-A A) grows faster 
than any polynomial as A —)■ +oo, and hence that the family {oj} must be infinite. 

We deduce that TTQy is the probability distribution of the independent sum 

Y = Y, Z, + Z, 

i>l 

where Zi ~ Exp(ai) for all i > 1 with 0 < ai <02 < ... < CLn < ■ ■ ■; and Z has a completely 
monotone density. This entails in passing that foy is real analytic on (0, + 00 ) - see also [7] for a 
weaker result in the framework of stochastic differential equations. Comparing the above infinite 
sum with the sum of n + 1 independent copies of Zn+i also shows that vroy(0, z) = O(z^) for every 
n > 1 as z —>■ 0+ (this latter fact is known in the literature as Ray’s estimate, see Problem 4.6.5. 
p.l34 in 0 ) and hence that /oy^( 0 +) = 0 for every n > 0 . 

The remainder of the proof is close to that of the Theorem in m- Introduce the notation 

_ f + if n is even 
I — if n is odd 

for every n G N. For {ei,..., e„,} some finite sequence in {—,0, +}, we say that a smooth function 
/ : ( 0 , + 00 ) —>■ M is of type ei... if it has limits (hnite or infinite) at zero and at inhnity, 
vanishes on a hnite set, and if the ordered sequence of its signs on [ 0 , + 00 ] is given by {ei,..., e^,}. 
For brevity, we will write / ~ ei... to express this property. With this notation, we need to 
show that 


( 2 ) 


An) 

JOy 




0 ±° 0 ±^ 0 --- ±” 0 , 


n > 1. 
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Observe first that this property holds true for n = 1, by Rosier’s result and the fact that f^y has 
isolated zeroes on (0,+oo). To show the property for all n > 2 we will use the same argument as 
in m- Fix n >2 and consider the independent decomposition 

y = + Yn, 

with 

y^ = ^ Z* + y and £ Y, 

i=l i>n-\-3 

Setting gn and hn for the respective smooth densities of Yn and 1^, we borrow from the Proposition 
in na the following key-result. 

Proposition f|ll]b With the above notations, one has 

(3) ~ 0±°0±^0---±*0, 0<i<n + l. 

In particular, the function has a extension on M and we can differentiate the convolution 
product to obtain 

poo 

9Y^\v)hn{u-v)dv, u>0. 

Using now the standard notation of [5], set S~{f) resp. S^{f) for the number of sign changes of a 
real function / on (0, -boo), the zero terms being discarded resp. included. By ([3|), Theorem 3.1.(a) 
p. 21 in [5] and the fact that the function is a Polya frequency function of infinite order - see 
e.g. Example 3.2.2 in [2], we have 

S-(/£*'') < S-(s<"+‘>) = n + l. 

This inequality entails first that fl^\+OQ) — 0. Indeed, it is easily seen from the fact that fiyy is a 

(n) 

density function that 0 must be a limit point of f^y' at infinity. If it is not a true limit, then Rolle’s 
theorem entails that S~{fQy~^^^) = -boo, a contradiction. 

Since (0-b) = (+oo) = 0 for all n > 0, an induction based on Rolle’s theorem shows first 

that S~^{fQy^) > n for all n > 1. On the other hand, since has isolated zeroes on (0,-boo) and 
vanishes at zero and infinity, again Rolle’s theorem shows that 

S+(/i^)) < S-{f^Y^)-l < n 

for all n > 1. This shows *S'^(/oy^) = n for all n > 1, which readily entails ([2]) as required. 

2.2. Proof of the only if part. Again, we can suppose x = 0, translating m if necessary. If m 
has exactly n points of increase on (0, y) then an induction based on ([1]) shows that A i-A- V'(y) A) is a 
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polynomial of degree n and hence, that TToy is the probability distribution of the hnite independent 
sum 

n 

Y = Y, Z, + Z, 

i=l 

with the above notation. It follows then from the Proposition in m that 

~ ±° 0 ±^ 0---±”0 

for all i > n, so that foy is not bell-shaped. 

2.3. Final remark. When m has exactly n points of increase on {x,y), the Proposition in |11] 
also provides some visual insight on the density fx,y When n = 1 for example, it shows that fx,y 
is ’’whale-shaped” (concave and then convex, all derivatives vanishing once). 
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